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Abstract 

We discuss an eigenvalue problem which arises in the studies of asymptotic stability of a 
self-similar attractor in the sigma model. This problem is rather unusual from the viewpoint 
of the spectral theory of linear operators and requires special methods to solve it. One of 
such methods based on continued fractions is presented in detail and applied to determine 
the eigenvalues. 

1 Introduction 

Many nonlinear evolution equations have the property that solutions which are initially smooth 
become singular after a finite time. The nature of this phenomenon, usually called blowup, 
has been a subject of intensive studies in many fields ranging from fluid dynamics to general 
relativity. The problem whether the blowup can occur, and if so, what is its character, is 
very difficult for some major evolution equations in physics, such as Navier-Stokes equations or 
Einstein's equations. Thus, in order to get some insight, it seems useful to study toy models. 
This paper is concerned with such a toy model, a nonlinear radial wave equation 

2 sin(2u) , , 

u t t - u rr - -u r H — — = 0, (1) 

where r is the radial variable and u = u(t,r). Equation Q describes equivariant wave maps 
from the 3 + 1 dimensional Minkowski spacetime into the three-sphere (see PP for the derivation). 
In the physics literature this model is usually referred to as the sigma model. 

The central question for equation Q is whether solutions starting from smooth initial data 

u(0,r)=f(r), u t (0,r)=g(r) (2) 

may become singular in future? A hint towards answering this question comes from the scaling 
argument. Note that equation JTJ is scale invariant: if u(t,r) is a solution, so is u\(t,r) — 
u(t/X, r/X). Under this scaling the conserved energy 

oc 

E[u] =J^u 2 t +u 2 r + 2 -^P^j r 2 dr (3) 
o 

transforms as a homogenous function of degree one: -B[ua] = XE[u], which means that equation 
(|T)) is supercritical in the language of the theory of nonlinear partial differential equation. For 
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supercritical equations it is energetically favorable that solutions shrink to small scales so singu- 
larities are expected to develop from sufficiently large (in a suitable norm) initial data. Although 
there are no rigorous results in this respect for equation (y , an explicit example of a singularity 
forming from smooth initial data is known. This example, first pointed out by Shatah [5] and 
later found in closed form by Turok and Spergel [3] is provided by the self-similar solution 

T 

u(t,r) = Uq(p) = 2arctan(p), where p = — — - (4) 
and T > is a constant 1 . Since 

drU (p) = (5) 

r— 1 — I 

the solution Uo(p) becomes singular at the center when t f T. By the finite speed of propagation, 
one can truncate this solution in space to get a smooth solution with compactly supported initial 
data which blows up in finite time. 

In fact, the self-similar solution Z7 is not only an explicit example of singularity formation, but 
numerical simulations indicate that it appears as an attractor in the dynamics of generic initial 
data pQ. We conjectured in [TJ that generically the asymptotic profile of blowup is universally 
given by Uq, that is 

lim u(t, (T — t)r) = Uo(r). (6) 

To prove this conjecture one needs to understand the mechanism responsible for the process of 
local convergence to the self-similar solution Uq. Such a mechanism is relatively well-understood 
for nonlinear diffusion equations where the global dissipation of energy is responsible for the 
convergence to an attractor, however very little is known for conservative wave equations where 
the local dissipation of energy is due to dispersion. 

In this paper, as the first step towards proving the above conjecture, we describe in more 
detail how the limit 10 is attained. To this end, in section 2 we consider the problem of linear 
stability of the solution Uq. This leads to an eigenvalue problem which is rather unusual from 
the standpoint of spectral theory of linear operators. In section 3 we solve this problem using 
the method of continued fractions. Finally, in section 4 we present the numerical evidence that 
the deviation of the dynamical solution from the self-similar attractor is asymptotically well 
described by the least damped eigenmode. 



2 Formulation of the eigenvalue problem 

In order to analyze the problem of linear stability of the self-similar solution Uq it is convenient 
to define the new time coordinate r = — ln(T — t) and rewrite equation in terms of U(t, p) — 
u(t, r) 

U TT + U T + 2p U pT - (1 - p 2 )(U pp + -U p ) + = 0. (7) 

In these variables the problem of finite time blowup in converted into the problem of asymptotic 
convergence for r — > oo towards the stationary solution Uq (p) . Following the standard procedure 
we seek solutions of equation Q in the form U(r,p) — Uo(p) + w(r,p). Neglecting the 0(w 2 ) 
terms we obtain a linear evolution equation for the perturbation w(t, p) 

2w 2 , 2cos(2J7 ) 

W TT + W T + 2/9 Wp T - (I - p ){W„„ + -W„) H 5 W = 0. (8) 

P P 

Uq is the ground state of a countable family of self-similar solutions U n (n = 0, 1, . . . ) 0. However, all n > 
solutions are unstable so they do not appear in the dynamics for generic initial data. 
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Substituting w(t, p) — e Xr v(p)/p into © we get the eigenvalue equation 

-(1 - p 2 )v" + 2Xpv' + A (A - l)v + ^-v = 0, (9) 

p 2 

where 

V{p) = 2 cos(4 arctanp) = 2 ^~ 6p - ■ (10) 

p z (l + /r) z 

We consider equation (J5J) on the interval < p < 1, which corresponds to the interior of the past 
light cone of the blowup point (t = T, r = 0). Since a solution of the initial value problem for 
equation starting from smooth initial data remains smooth for all times t < T, we demand 
the solution v(p) to be analytic at the both endpoints p = (the center) and p = 1 (the past 
light cone). Such a globally analytic solution of the singular boundary value problem can exist 
only for discrete values of the parameter A, hereafter called eigenvalues. 

A straightforward way to find the eigenvalues would be to use the Frobenius method. The 
indicial exponents at the regular singular point p = are 2 and —1, hence the solution which is 
analytic at p = has the power series representation 

oo 

Since there are no complex singularities in the open disk of radius 1 about p = 0, the series (|11|) 
is absolutely convergent for < p < 1. At the second regular singular point, p = 1, the indicial 
exponents are and 1 — A so, as long as A is not an integer (below we shall discuss this case 
separately), the two linearly independent solutions have the power series representations 

oo oo 

«i(p) = E ti 1 ^ 1 - p) n i v ^p) = E ^ 2) a - p) n+1 ~ x - (12) 

These series are absolutely convergent for < p < 1. If A is not an integer then only the solution 
v\(p) is analytic at p = 1. From the theory of linear ordinary differential equations we know 
that the three solutions Vo(p), vi(p), and Wa(p) are connected on the interval < p < 1 by the 
linear relation 

v (p) = A(X) Vl (p) + B(\)v 2 (p). (13) 

The requirement that the solution which is analytic at p = is also analytic at p = 1 serves as 
the quantization condition for the eigenvalues B(X) = 0. Unfortunately, the explicit expressions 
for the connection coefficients ^4(A) and B(X) are not known for equations with more than three 
regular singular points. In the next section we shall present an indirect method which goes 
around this difficulty. 

We remark in passing that alternatively the eigenvalues can be computed numerically using 
a shooting-to-a-midpoint technique. With this technique one approximates the solutions Vq(p) 
and Vi(p) by the power series (|ll|l and (|12fl . truncated at some sufficiently large n, and then 
computes the Wronskian of these solutions at a midpoint, p = 1/2, say. The zeros of the 
Wronskian correspond to the eigenvalues. Although this technique generates the eigenvalues 
with reasonable accuracy, it is computationally very costly, especially for large negative values 
of A, because the power series ()ll|l and (|12f) converge very slowly. Note also that shooting 
towards p = 1 fails completely for large negative A because the solution i>2 (p) is subdominant at 
p = 1, that is, it is negligible with respect to the analytic solution v\(p). 
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3 Solution of the eigenvalue problem 



In this section we shall solve the eigenvalue problem @ using a method which exploits an inti- 
mate relationship between recurrence relations and continued fractions. Although this method 
is not widely known, it is in fact quite old and has been applied in the past to determine the 
bound states of the hydrogen molecule ion |5] and quasinormal modes of black holes jo]. 

The key idea is to determine the analyticity properties of the power series solution Vo(p) from 
the asymptotic behavior of the expansion coefficients a n . In order to implement this idea it is 
convenient to change the variables 

v(p) = (2-x)^y(x), x = ^^. (14) 

In terms of these variables equation ljUJ) takes the form 

x 2 {\ - x){2 - x)y" + x[l - (1 + X)x(2 - x)]y' - ^[X 2 x(l - x) + 9x 2 - 17x + 4}y = 0. (15) 

The reason of making the transformation (|14Jl is twofold. First, the transformation of the 
independent variable rearranges the singular points of equation 10 in such a way that, without 
changing the points p = and p = 1, moves the bothersome singularities at p = H (which lie 
on the unit disk around p = and obstruct the analysis of analyticity of the power series (|11|) 
at p = 1) to infinity and moves p — oo to x — 2. Second, by factoring out the singular behavior 
at x = 2 the number of terms in the recurrence relation for the coefficients of the power series 
solution around x = is reduced from four to three. The indicial exponents at x — are 1 and 
— 1/2 so the solution which is analytic at x — has a power series expansion 

oo 

y (x) = a « a; "' a i^°- ( 16 ) 

n=l 

Substituting this series into equation l|15[) we get the three-term recurrence relation 
p 2 (0)a 2 +pi(0) ai = 0, 

P2(n)a n+2 +pi(n)a n+1 + p Q (n)a n = 0, n = l,2, ... (17) 

with the initial conditions eto = and a\ = 1 (normalization) where 

p 2 (n) = 8n 2 + 28n + 20, (18) 
pi(n) = -12n 2 -(20 + 8A)n-A 2 -8A + 9, (19) 
p (n) = An 2 + 4An + A 2 - 9. (20) 

The series (jlfijl is absolutely convergent for < x < 1 and in general is divergent for x > 1. In 
order to determine the analyticity properties of the solution yo(x) at x = 1 we need to find the 
large n behavior of the expansion coefficients a n . The three-term recurrence relation (|17|l can 
be viewed as the second order difference equation so it has two linearly independent asymptotic 
solutions for n — » oo. Following standard methods (see, for example, 0) we find 

oo (1) oo (2) 

^W^E^' and «[^2-Vl^^, (21) 

where the coefficients Cg can be determined recursively from l|17|) . Thus, in general the solution 
of the recurrence relation (|17fl behaves asymptotically as 

a„~ Cl (A)a«+c 2 (A)4 2) . (22) 



4 



If the coefficient ci(A) is nonzero then 



dn+l a n+l 1 / 00 \ 

1 1 as n — > oo, (23) 



hence the power series (1 1 01) is divergent for cc > I (in fact, it has a branch point singularity at 
x — 1). On the other hand, if ci(A) = then 



(2) , 

a-n+l a n+l _j_ 



as n — > oo, (24) 



and the power series 1(16(1 is absolutely convergent for x < 2, in particular the solution yoC^) is 
analytic at x = 1. 

The advantage of replacing the quantization condition B(X) — in the connection formula 
(113(1 by the equivalent condition ci(A) = follows from the fact that ci(A) is the coefficient of 
the dominant solution in 1)22(1 . in contrast to B(X) which is the coefficient of the subdominant 
solution in 113(1 . In the theory of difference equations (recurrence relations) a subdominant 
solution, that is a solution which is asymptotically negligible with respect to any other solution, 
is called the minimal solution. In contrast to a dominant solution, the minimal solution, if it 
exists, is unique. The condition ci(A) = is thus equivalent to the requirement that the solution 
of the recurrence relation 1(17(1 starting with ao = and a\ = 1 is minimal. In order to find 
when this minimal solution exists we shall use now a relationship between three-term recurrence 
relations and continued fractions. Let 



A n = P -^, B n = P 4±, r n =^. (25) 
...... „ 



Then, we can rewrite ((1711 as 



Pa(n)' P2(n) 



A n + r n+1 

and applying this formula repeatedly we get the continued fraction representation of r n 

B n B n+ i B n+ 2 
r n = ... 

A n — A n +\— A n +2 — 



(26) 



(27) 



Pincherle's theorem |7] says that the continued fraction on the right hand side of equation ((27(1 
converges if and only if the recurrence relation 117(1 has a minimal solution a™ m and, moreover, 
in the case of convergence, equation 1(27(1 holds with r n = a"^ r l/a™ ln for each n. 

Using Pincherle's theorem and setting n — 1 in l(27|) we obtain the eigenvalue equation 

The continued fraction in 1(28(1 . which by Pincherle's theorem is convergent for any A, can be 
approximated with essentially arbitrary accuracy by downward recursion starting from a suf- 
ficiently large n — N and some (arbitrary) initial value rjy. The roots of the transcendental 
equation l(2"%|l are then found numerically (see table 1). 

A glance at table 1 shows some interesting properties of the spectrum. First, all the eigen- 
values are real. Second, there are no eigenvalues A > 1. Third, there are two integer eigenvalues 
A = 1 and A = — 2. Below we discuss these properties in detail. 
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n 





1 


2 


3 


4 


5 


An 


1 


-0.542466 


-2 


- 3.398382 


-4.765079 


-6.102295 


n 


6 


7 


8 


9 


10 


11 


An 


-7.297807 


-7.765347 


-8.853889 


-10.1228208 


-11.196495 


-11.802614 



Table 1: The first twelve eigenvalues. 



Why the eigenvalues are real? 

We find this property surprising because we cannot see any a priori reason which forbids complex 
eigenvalues. On the contrary, it is easy to deform the potential V in equation JJjJ, without chang- 
ing the character of singularities at the endpoints, in such a way that the resulting eigenvalues 
are complex. Thus, the reality of eigenvalues is related to the special form of the potential V. 
It remains a puzzle whether this relationship is accidental or there is something deep to it. 



The eigenvalue A = 1 

This eigenvalue is due to the freedom of changing the blowup time T. To see this, consider a 
solution Uo( T J~_ t ) with shifted blowup time. In terms of the similarity variables r = — ln(T — t) 
and p = r / (T — t) we have 

u ° (r^t) = Uo (irbO ■ where e = T '- T - (29) 

Thus, the perturbation induced by the shift of blowup time has the form 

w{p,r) = -ee T U' Q {p), (30) 

and consequently the mode 

«(p) = PK(P) = T~7~~2 ( 31 ) 
1 + p z 

solves equation with A = 1. Since this mode is evidently analytic at both p = and p = 1, 
we conclude that A = 1 is an eigenvalue. We emphasize that this positive eigenvalue should not 
be interpreted as the physical instability of the solution Uq, as it is an artifact of introducing 
the similarity variables and does not show up in the dynamics for u(t,r). 



Nonexistence of eigenvalues A > 1 

Using the transformation 

v(p) = (l-p 2 )-$z(p), (32) 
we can put equation into a standard Sturm-Liouville form 



cP_ 1-p 2 
dp 2 + p 2 



Az = f iz, where A = -(1 - p 2 ) 2 ^ + ^~j-V(p), fi = A(2 - A). (33) 
In the space of functions 



V = L 2 ([0,1], (// ' 



(i- P 2 y 
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the operator A is self- adjoint so the eigenvalues p are real. Both endpoints are of the limit-point 
type. Near p = an admissible (that is, belonging to V) solution behaves as z(p) ~ p. Near 
p = 1 two independent solutions behave as 

Mp)~(l"P)* (1±VI=Ji) , (34) 

so only the solution z + with /i < 1 is admissible. Now, we shall show the operator A has no 
eigenvalues. To see this, note that the solution with fi = 1, corresponding to the gauge mode 



z(p) = (35) 
1 + p- 

has no zeros. This implies by a standard theorem for Sturm-Liouville operators that there are 
no eigenvalues below p = 1. Since p = 1 is not an eigenvalue (because the mode (|35|l does not 
belong to T>\ we conclude that the operator A has the purely continuous spectrum p > 1. 

What does this fact tell us about the eigenvalues of our problem? Rewriting Ij34[) in terms 
of A we have 

Mp)~(l-p)*, MpWl-p) 1 -^, (36) 

so comparing with (|32[1 we see that only z + leads to a solution v(p) which is analytic at p = 1. 
For Re(X) > 1 the solution z + belongs to 2? so in this case there is one-to-one relationship 
between the eigenvalues of the operator A and the eigenvalues of our problem. Since the former 
has no eigenvalues, we infer that our problem has no eigenvalues with Re(X) > 1. This result 
was obtained previously in 0]. 

We point out that for Re(X) < 1 the requirements of square-integrability and analyticity 
near p = 1 are mutually exclusive so in this case there is no relationship between between the 
eigenvalues of the operator A and the eigenvalues of our problem. 

The algebraically special eigenvalue A = —2 

If 1 — A = TV is a positive integer then the two linearly independent power series solutions of 
equation l|15H around x = 1 are 



Vi 



{x) = Y^b£\l-xY+ N , y 2 {x)=C N y 1 {x)Hl-x) + Y d bW{l-x) n . (37) 



The solution y\{x), corresponding to the larger indicial exponent, is analytic at x = 1 but the 
solution 2/2(a0, corresponding to the smaller indicial exponent, involves a logarithm in general. 
However, for TV = 3 we have an exceptional case: the coefficient C3 vanishes and both solutions 

are analytic at x = 1. The best way to see this is to solve the recurrence relation for the 

(2) 

coefficients 6„ assuming temporarily that the solution j/2 contains no logarithm. We find 

5(2) = P2n-l(\) , s 

A(A + l)(A + 3)...(A + n-l)' [ ' 

where P2n-i(A) is a polynomial in A of order In — 1. Moreover this polynomial has no integer 

roots. When 1 — A = N is a positive integer different from 3, then the expansion coefficient 
(2) 

b N does not exist which contradicts the assumption that j/2 contains no logarithm. The only 
exception is TV = 3 because the denominator in (|38|l has no (A + 2) term. Thus, for A = — 2 
the singularity at x = 1 is apparent - each solution which is analytic at x = is automatically 
analytic at x — 1 as well. This proves that A = —2 is the eigenvalue. 



7 



4 Numerical verification 



According to the linear stability analysis presented above the convergence of the solution u(t, r) 
towards the self-similar attractor Uq should be described by the formula 



u(t,r) = U(r,p) = U (p) + } c k e x « T v k (p)/p ~ U (p) + c 1 e MT v l (p)/p 



k=l 



as 



oo, 



(39) 



where Vk(p)/p are the eigenmodes corresponding to the eigenvalues and Ck are the expansion 
coefficients. In order to verify (|39[1 we solved equation numerically for large initial data 
leading to blowup, expressed the result in the similarity variables, and computed the deviation 
from Uq for t f T. The result (see figure 1) shows that, in perfect agreement with the formula 
J2HJ), the deviation from C/ is described by the least damped eigenmode v\. This makes us feel 
confident that the calculation presented in section 3 contains no algebraic errors. 




1.2 1.4 



Figure 1: We plot the deviation of the dynamical solution u(t,r) from the self-similar solution Uo = 
2 arctan(p) at some moment of time close to the blowup time. The solid line shows the least damped 
eigenmode ci(T — t)~ Xl vi(p) / p, where the coefficient ci is fitted once for all times. 

We remark that the formula (|39|l could be used to compute the eigenvalues numerically directly 
from the dynamics rather than by solving the eigenvalue equation. Such a computation was 
performed by Donninger [S] with the result which is in rough agreement with table 1 (rough, 
as the dynamical computation of eigenvalues is by far less accurate than the continued fraction 
method). 
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